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G\ : Abstract 

f-^. ■ The spectrum of the Hermitian Hamiltonian ^p 2 + ^m 2 x 2 + gx 4 (g > 0), 

in 

which describes the quantum anharmonic oscillator, is real and positive. The 

00 

non-Hermitian quantum- mechanical Hamiltonian H = ^p 2 + |m 2 x 2 — gx 4 , 



where the coupling constant g is real and positive, is 'PT-symmetric. As a 
consequence, the spectrum of H is known to be real and positive as well. Here, 



Oh 

• it is shown that there is a significant difference between these two theories: 



When g is sufficiently small, the latter Hamiltonian exhibits a two-particle 



bound state while the former does not. The bound state persists in the 
corresponding non-Hermitian PT-symmetric — gcf) 4 quantum field theory for 
all dimensions < D < 3 but is not present in the conventional Hermitian 
gcj) A field theory. 

In this Letter we show that the spectrum of the non-Hermitian[| PT-symmetric quartic 
Hamiltonian 



1 Contrary to appearances, this Hamiltonian is not Hermitian because its eigenfunctions are re- 
quired to obey boundary conditions in the complex plane [see (||)]. 



H = ^p 2 + ^m 2 x 2 -gx 4 {g > 0) (1) 

possesses bound states for sufficiently small g. Moreover, when the Euclidean space-time 
dimension D satisfies < D < 3, the corresponding PT-symmetric quantum field theory 
described by the Lagrangian 

C = \{V<PY+ l -m 2 <f-g^ (g>0) (2) 

also possesses bound states for small g. The related conventional Hermitian gx 4 and g<f) 4 
theories with g > do not possess such bound states. 

Quantum-mechanical Hamiltonians that are non-Hermitian but possess VT symmetry 
have recently been studied in great detail by many authors PHI^fl. It is known that the 
entire spectrum of the PT-symmetric anharmonic oscillator is real and positive. The first 
proof of the reality of the spectrum is given in Ref. ||13|| . Direct numerical evidence for the 



reality and positivity of the spectrum can be found by performing a Runge-Kutta integration 
of the associated Schrodinger equation [E. Alternatively, the large-energy eigenvalues of 
the spectrum can be calculated with great accuracy by using conventional WKB techniques 
|l|j3|JT5|] . An easy way to demonstrate the reality and positivity of the spectrum is to calculate 
exactly the spectral zeta function (the sum of the inverses of the eigenvalues). This was done 
by Mezincescu and Bender and Wang |7||§. 

The calculations that are performed in this paper are perturbative; that is, quantities 
are determined as formal power series in the coupling constant g. Some explanation must 
be given to justify the application of perturbative methods to PT-symmetric theories. We 
will argue that to obtain a perturbation expansion for the theories described in ([]]) and 
(0) we need only calculate the perturbation series for the corresponding quantity in the 
conventional Hermitian gx 4 and g<j) 4 theories and change the sign of g in these expansions. 
This procedure is correct but it is nontrivial, as we now argue. 

First, consider the quantum- mechanical theory in (|l|). The corresponding differential 
equation for the energy eigenfunctions ip(x) is 
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— ~if>"(x) + -m 2 x 2 ip(x) — gx 4 ijj(x) = Eip(x) (g > 0), (3) 

where the wave function ip(x) satisfies boundary conditions in two wedges of 60° opening 
angle in the complex- x plane [I]: 

lim ip(x) = (— 7r/3 < argx < and — n < argx < — 27r/3) . (4) 

\x\— >oo 

These two wedges lie just below the real axis. Because they are symmetric with respect to 
the imaginary axis they enforce the VT symmetry of the problem. In the limit as g — > 0, 
the differential equation in (0) becomes that of the harmonic oscillator and the boundary 
conditions are imposed in two wedges of 90° opening angle in the complex-x plane: 

lim if)(x) = (— tt/4 < argx < 7r/4 and — 57r/4 < argx < — -tt) . (5) 

Because the wedges for g > and the wedges for g = have a region in common, conventional 
perturbative methods are justified. 

However, one must be extremely careful to avoid drawing wrong conclusions on the basis 
of perturbative calculations. For example, a perturbative calculation of the one-point Green's 
function (the expectation value of x in the ground state) gives to all orders in powers of 
the coupling constant g. Nevertheless, the correct result is nonzero, pure imaginary, and 
exponentially small ||16|| : 



^ f dxXlbZ(x) . rJ , 0/9 n/o -4m 3 

where K = 8^3'^ /T 2 (1/3) is a positive constant. The nonperturbative result in (||) is 
due to a soliton solution of the classical field equations. 

There are even more subtle effects in which one may think that there are nonpertur- 
bative effects when there are none. For example, consider the cubic anharmonic oscillator 
Hamiltonian 

H = -p 2 + -m 2 x 2 + gx 3 (g real). (7) 
A perturbative calculation of the ground-state energy has the form 



1 oo 

E ~-m + Y,an(g 2 T fo-O). (8) 
A i 

The coefficients a n in this series are real and all have the same sign and thus the series is 
not Borel summable. The nonsummability of this series indicates the presence of a branch 
cut in the complex-^ plane on the real-g axis. Thus, the function E (g) is complex when g is 
real. The imaginary part of Eo(g) is a nonperturbative effect; it is exponentially small when 
g is real and small. The energies of the cubic oscillator are complex because of tunneling 
effects on the real axis. 

The corresponding PT-symmetric cubic oscillator, whose Hamiltonian is 

H = -p 2 + -m 2 x 2 + igx 3 ((/real), (9) 

is obtained from the Hamiltonian in (0) by replacing g with ig. This replacement makes the 
series in (|J) Borel summable and verifies the result that eigenvalues of H in (|)|) are 



all real. 

However, this argument is quite tricky. Let us apply the same reasoning to the conven- 
tional Hermitian quantum-mechanical anharmonic oscillator described by 

H = ^p 2 + ~m 2 x 2 + gx* (g > 0). (10) 

The Rayleigh-Schrodinger perturbation series for the ground-state energy of this Hamilto- 
nian has the form 

oo 

Eo^-m + J £an9 n G?^0+), (11) 
A i 

where the coefficients a n are all real and alternate in sign. Hence, this series is Borel 
summable and we conclude (correctly) that the ground-state energy is real when g > 0. 
We correctly obtain the perturbation expansion for the ground-state energy of the non- 
Hermitian PT-symmetric anharmonic oscillator ([!]) by replacing g with —g. Now, the 
coefficients in perturbation series for the ground-state energy no longer alternate in sign; 
the coefficients all have the same sign (they are all negative). However, all of the eigenvalues 



of the Hamiltonian ([J) are real. Even though the series is not Borel summable, the imaginary 
part of the ground-state energy is exactly zero due to the presence of the soliton that leads 
to the nonperturbative result in (^). 

The same discussion applies to the non-Hermitian PT-symmetric quantum field theory 
in (|2|). The perturbation expansions of the Green's functions of this theory are obtained 
using the same Feynman rules as for the conventional g0 4 quantum field theory, with the 
one change that the vertex amplitude —24g is replaced with 24g. However, again one must 
be very careful about nonperturbative effects. For example, the one-point Green's function, 
which vanishes in the conventional g<p A quantum field theory, is now nonvanishing, pure 
imaginary, and exponentially small in the — g0 4 quantum field theory [16[ . 



Having made these preliminary remarks, we turn to the question of bound states in VT- 
symmetric theories. Let us first examine a two-particle state for the conventional quantum 
anharmonic oscillator Hamiltonian (|T0|). The small-g Rayleigh-Schrodinger perturbation 
series for the fcth energy level of this Hamiltonian is ]17| 
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k + ^ + ^(2k 2 + 2k + l)e + 0(e 2 ) 



0+), (12) 



where the dimensionless expansion parameter is e = g/m 3 . In this theory the renormalized 
mass M is defined as the first excitation above the ground state. Thus, from (|12|), 



M = E 1 - E ~ m [l + 3e + 0(e 2 )] (e + ). (13) 

To determine if the two-particle state is bound, we examine the second excitation above 
the ground state using (0). We define 

B 2 = E 2 -E ~m [2 + 9e + 0(e 2 )] (e -> + ). (14) 

If B 2 < 2M, then a two-particle bound state exists and the binding energy, which is a 
negative quantity, is B 2 — 2M. If B 2 > 2M, then the second excitation above the vacuum 
is interpreted as an unbound two-particle state. We can see from ( |T4"D that in the small- 
coupling region, where perturbation theory is valid, the conventional anharmonic oscillator 



does not possess a bound state. Indeed, using a variety of methods (WKB theory, variational 
methods, numerical calculations) one can show that there is no two-particle bound state for 
any value of g > 0. Because there is no bound state the gx 4 interaction may be considered 
to represent a repulsive force.Q 

We obtain the perturbation series for the PT-symmetric oscillator, whose Hamiltonian 
H is given in (|l|), from the perturbation series for the conventional anharmonic oscillator by 
replacing e with — e. Thus, while the conventional anharmonic oscillator does not possess a 
two-particle bound state, evidently the PT-symmetric oscillator does indeed possess such a 
state. The binding energy of this state can be measured in units of the renormalized mass 
M. For this purpose we define the dimensionless binding energy A 2 by 

A 2 ^ 2 ^ 2M ~-3e + 0(6 2 ) ( e -0 + ). (15) 

In Fig. [TJ, A 2 is plotted as a function of the dimensionless coupling constant e. Observe 
that the bound state disappears when the value of e increases beyond e = 0.0465 .... As e 
continues to increase, A 2 reaches a maximum value of 0.427 at e = 0.13 and then approaches 
the limiting value 0.28 as e — > oo. 

It is interesting that the bound state for the PT-symmetric anharmonic oscillator occurs 
only for values of e less than about 0.0465. The following heuristic argument helps to explain 
why the characteristic size of e is so small. For the Hamiltonian in (0) the potential has 
the form V(x) = |m 2 x 2 — gx 4 . On the real-x axis, the maximum value of V(x) occurs at 
x 2 = \m 2 /g and at this point V mabX = T6 m V<?- In order to have a two-particle bound state 

2 In general, a repulsive force in a quantum field theory is represented by an energy dependence in 
which the energy of a two-particle state decreases with separation. The conventional anharmonic 
oscillator Hamiltonian corresponds to a field theory in one space-time dimension, where there 
cannot be any spatial dependence. In this case the repulsive nature of the force is understood to 
mean that the energy 5 2 needed to create two particles at a given time is more than twice the 
energy M needed to create one particle. 
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it is necessary for the height of the potential to be high enough to bind three states, the 
largest of which is of order |m. Setting this value equal to V m3iX gives e = which is near 
the value for which the two-particle state is most strongly bound. Note that the potential 
V(x) does not allow states to tunnel out of the well, even though the potential decreases for 
large real x. This is because the PT-symmetric boundary conditions are not imposed on 
the real axis but rather in a pair of wedges in the complex-x plane that lie below the real 
axis [see Eq. (|J)]. 

In the PT-symmetric anharmonic oscillator, there are not only two-particle bound states 
for small coupling constant but also fc-particle bound states for all k > 2. To see why, we 
calculate the energy of the kth excitation above the vacuum: 



B k = E k - E ~ m 



k- l(k 2 + k)e + 0(e 2 ) 



0+), (16) 



which is the generalization of (|T4|). Thus, the dimensionless binding energy is 



A fc = ^^~-^(A;-l) e + 0( e 2 ) ( e -0+), (17) 

This equation reduces to flT5| ) when k = 2. The key feature of this equation is that the 
coefficient of e is negative. Since the dimensionless binding energy becomes negative as e 
increases from 0, there is a fc-particle bound state. The dimensionless binding energies for 
the first five A;-particle bound states are shown in Fig. |2|. Note that the higher multiparticle 
bound states cease to be bound for smaller values of e; starting with the three-particle 
bound state, the binding energy of these states becomes positive as e increases past 0.039, 
0.034, 0.030, and 0.027. Furthermore, the heuristic argument given above for the value of e 
at which the multiparticle states are most strongly bound is qualitatively correct; one can 
see that the fc-particle bound state is most strongly bound at a value of e that is roughly 
1/(16*; + 8). 

Figure |2] shows that for any value of e = g/m 3 there are always a finite number of bound 
states and an infinite number of unbound states. The number of bound states decreases 
with increasing e until, as we see in Fig. [I], there are no bound states at all. Observe that 



there is a range of e for which there are only two- and three-particle bound states. (This 
situation is analogous to the physical world in which one observes only states of two and 
three bound quarks.) In this range of e if one has an initial state containing a number 
of particles (renormalized masses), these particles will clump together into bound states, 
releasing energy in the process. Depending on the value of e, the final state will consist 
either of two- or of three-particle bound states, whichever is energetically favored. Note also 
that there is a special value of e for which two- and three-particle bound states can exist in 
thermodynamic equilibrium. 

Let us now generalize the result for the quantum-mechanical two-particle bound state 
in fll5|) to the case of the D-dimensional PT-symmetric quantum field theory in (|2|). To 
show that there is a bound state, we need only demonstrate that to leading order in the 
coupling constant, the binding energy becomes negative as the coupling constant increases 
from 0. The Feynman rules for the calculation are the conventional ones, except that the 
sign of the vertex is reversed: 24g for a four-point vertex amplitude and l/(p 2 + M 2 ) for a 
line amplitude. Note that to leading order in the coupling constant we use the renormalized 
mass M rather than the unrenormalized mass m in the propagator amplitude and we exclude 
self-energy loops in the graphs. 

To calculate the energy of the bound state to leading order we sum up all "sausage-link" 
graphs and identify the bound-state pole in the geometric series. These graphs are shown in 
Fig. |3|. The result for the dimensionless binding energy to leading order in the dimensionless 
coupling constant e = gM D ~ A is 



where D is the Euclidean space dimension. The formula reduces to that in for the 
case D = 1. This formula is valid for < D < 3. Note that we have performed a mass 
renormalization to leading order, but we have not done a wave function or coupling-constant 
renormalization in this calculation. 

We conclude by remarking that the behavior of a gcj) 3 theory is the reverse of that of a 




(18) 
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g0 4 theory. A g<f) 3 theory represents an attractive force. The bound states that arise as a 
consequence of this force can be found by using the Bethe-Salpeter equation. However, the 
g<p 3 field theory is unacceptable because the spectrum is not bounded below. If we replace 
g by ig, the spectrum now becomes real and positive, but the force becomes repulsive and 
there are no bound states. The same is true for a two-scalar theory with interaction of the 
form ig(f> 2 X- This latter theory is an acceptable model of scalar electrodynamics, but has no 
analog of positronium. 
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FIGURES 
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FIG. 1. Dimensionless binding energy A2 of the two-particle bound state of the PT-symmetric 
anharmonic oscillator, whose Hamiltonian H is given in Eq. (1), plotted as a function of the 
dimensionless coupling constant e = g/m 3 . The quantity A2 represents the binding energy B2 
measured in units of the renormalized mass [see Eq. (15)]. Note that for small e the slope of the 
curve is —3, which verifies the asymptotic result in (15). Observe that the bound state disappears 

when e increases beyond 0.0465 As e continues to increase, A 2 reaches a maximum value of 

0.427 at e = 0.13. Then A2 decreases and approaches the limiting value 0.28 as e — > 00. 
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FIGURE 2 
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FIG. 2. Dimensionless binding energies A2, A3, A4, A5, and Aq for the two-, three-, four-, 
five-, and six-particle bound states plotted as functions of the dimensionless coupling constant e. 
Note that the multiparticle bound states cease to be bound as e increases past 0.039, 0.034, 0.030, 
and 0.027. 
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FIGURE 3 
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FIG. 3. So called "sausage-link" graphs contributing to the leading-order value of the energy 
of the two-particle bound state for the 'PT-symmetric — gcj) 4 quantum field theory whose Euclidean 
Lagrangian is given in Eq. (2). 
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